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^ ■ Abstract. We present a new library of synthetic spectra based on the stellar atmo- 

^ I sphere code PHOENIX. It covers the wavelength range from 500 A to 55 000 A with a 

resolution of R=500 000 in the optical and near IR, R= 1 00 000 in the IR and A^ =0. 1 A 
in the UV. The parameter space covers 2 300 K < T^s < 8 000 K, 0.0 < log(g) < +6.0, 
-4.0 < [Fe/H] < +1.0 and -0.3 < [a/Fe] < +0.8. The librai-y is work-in-progress 
and going to be extended to at least T^ff = 25 000 K. We use a new self-consistent way 

tH- ■ of describing the microturbulence for our model atmospheres. The entire library of 

CK ! synthetic spectra will be available for download. 

Futhermore we present a method for fitting spectra, especially designed to work with 

ff^ [ the new 2nd generation VLT instrument MUSE. We show that we can determine stellar 

^D ■ parameters ( T^g, log(g), [Fe/H] and [a/Fe]) and even single element abundances. 



Keywords : spectral library 



H ■ 1. Introduction 



Presumably next year the Multi Unit Spectroscopic Explorer MUSE (iBacon et al.ll2010l) will see 



its first light as a second-generation instrument for ESO's Very Large Telescope at Paranal, Chile. 
The instrument is an highly efficient AO-supported integral field unit (IFU) and its outstanding 
combination of a large field of view (1x1 arcmin^) and high spatial sampling of 0.2" with a 
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spectral resolution of R=2000-4000 over the optical range of 4650-9300 A will allow us to obtain 
unprecedented observations. 

In our group, we intent to use it for the analysis of galactic globular clusters, which due to the 
heavy crowding towards the center are only accessible through their giants by other instruments. 
With MUSE, we will be able to point directly in the center of the clusters and obtain thousands 
of stellar spectra even from stars well below the main-sequence turnoff point with one single 
exposure. 

For the analysis of those spectra, we need a grid of model spectra that matches both the 
wavelength range and resolution of MUSE as well as our requirements for an extensive parameter 
space (as given by previous observations of globular clusters) and being able to adjust it as needed. 
Therefore a deci sion was made to create a n ew grid of model atmospheres and synthetic spectra 
with PHOENIX (JHauschildt & Baronll 19991) . 



We will present this new library together with a short description of the methods used for ana 



lyzing M USE spectra and some preliminary results on a simulated data cube. A paper (IHusser et al 
lin prep.) about our new synthetic stellar library with more detailed descriptions and informations 
for downloading the spectra is in preparation and will be published soon. 

2. Globular Clusters 

There are about 150 known globular clusters in our galaxy with masses of 10^ - 10^ Msun, which 
consist of very old stellar populations with an age of > lOGyr. A couple of years ago those 
populations were assumed to be very simple with a single isochrone. 

Unexpectedly, latest observation s showed evidence for the existence of multiple main-se- 



quences within globular clusters, e. g.pedin et al.l (120041) for Omega Centauri, and also for a split 
in the (sub)giant branch (iLee et al.l 19991) . This seems to be caused by multiple stellar populations 



with different abundances of helium and iron, at least for massive clusters. A variation of some 
lighter elements seems to be more ubiquitos, like the Na-O anti-correlation observed by Carretta 
(2009). Those observations indicate a complex enrichment history with multiple epochs of star 
formation. 

In globular clusters, we observe a velocity dispersion of 5 - 20 km/s and high central stellar 
densities of ~10^Msun/pc^- Given these high stellar densities, scen arios have been proposed for 



the fo rmation of intermediate-mass black holes in the cluster centres (IPortegies Zwart & McMillan 



20021) . An extrapolation of the tight relation between total mass and the mass of the black hole 
in the bulges of galaxies (Marconi and Hunt 2003) towards typical masses of globular clusters 
predicts black holes with a mass of ~10"' Msun- 

In contrast to field stars, where we observe a binary fraction of a bout 50%, in globula r clusters 
we find a value of only 30% or even less. Monte Carlo simulations (Ivanova et al. 2005) showed 
that the number of binaries in the core decreases rapidly over time. .Spitzer.(.1987.) already showed 
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Table 1. Parameter space of the grid. An extension in T^g up to 12 000 K is work in progress and up to 
25 000 K in planning. Alpha element abundances [a/Fe] i= are only available for 3 500 K < T^g < 8 000 K 
and -2 < [Fe/H] < 0. 



Range Step size 



T^eff [K] 


2 300- 


-7 000 


100 




7 000- 


-8 000 


200 


log(6f) 


0.0- 


--H6.0 


0.5 


[Fe/H] 


-4.0- 


--2.0 


1.0 




-2.0- 


--Hl.O 


0.5 


[a/Fe] 


-0.3- 


--H0.8 


0.1 



Table 2. Spectral resolution of the grid 


Rang 


e[A] 


Resolution 


500- 

3 000- 

25 000- 


-3 000 
-25 000 
-55 000 


AA 
R 
R 


= 0.1A 
« 500000 
« 100000 



that a depletion of binaries in the core is necessary for it to collapse, so core-collapsed clusters 
like M15 seem to be dynamically more evolved than others. A core-collapse could also result in 
the formation of new binaries. Therefore, studying the binary fraction in globular clusters is an 
important task in order to understand their evolution. 



3. The new PHOENIX grid 

The PHOENIX version 16 that we are using for calculating the grid uses a new equation of state 
called ACES, which is a state-of-the-art treatment of the chemical equilibrium in each layer of 
a stellar atmo s phere . The element abundances we used for the atmospheres were taken from 
Asplund et all i2009l) . 



The parameter space of the new PHOENIX grid that we are presenting in this paper is given 
in Table [1] An extension towards hotter stars including NLTE treatment of important elements is 
both work-in-progress (up to 12 000 K) and intended (up to 25 000 K). The grid is complete in its 
first three dimensions effective temperature T^-s, surface gravity log(gr) and metallicity [Fe/H]. 
Different alpha element (including O, Ne, Mg, Si, S, Ar, Ca and Ti) abundances are provided for 
3 500 K < Teff < 8 000 K and -2 < [Fe/H] < only. 

Although the primary intent for creating the grid was the analysis of MUSE spectra, we 
decided to increase both the wavelength range and the resolution (see Table |2]i, so that teams 
working on other existing and upcoming instruments will be able to use t hem for their purpo ses. 
Due to to th is, our spectra are applicable for the analysis of e. g. CRIRES (IKaeufl et alJ2004h and 
X-Shooter (iVemetet al.lboilb data. 



T.-O. Husser et al. 



stellar Mass [U^^^ 



















































































■■^■■■^^BB 






JJ^BBB^^Bj 


i"n 


■■^■■■^^■1 






■■^^■■^^■1 


HUU 



iiQDnnnnnnnnnnn 
on n D n n D n 



■ B B 

■ ^ □ 

EDO 

n n n 

H n D 

n n n 

n n n 

D n D 



Figure 1. Distribution of stellar masses for that part of the grid with solar abundances for different effective 
temperatures T^^ and surface gavities log(^). Color-coded is the stellar mass in units of solar mass from 
OMq (black) to 9Mq (white). 

In order to define a spherical symmetric atmosphere as it is used in PHOENIX, we need 
to define an effective temperature Teff, a surface gravity log(^) and either a radius ro or a mass 
Mi,. We decided to use the mass by taking a mass-luminosity relation L^^/Lq - {Mi^jMof for 
main-sequence stars and letting it tend towards higher values for giants and super giants: 



M^^c- M„ 



\5770K/ 



with values for the coefficient as given in the following table: 



(1) 



log(g) 



>4 >3 >2 >1.6 >0.9 >0 <0 

1 1.2 1.4 2 3 4 5 



Figure[T]shows the distribution of masses in our grid for solar abundances. 



PHOENIX uses the mixing length theory (IPrandtll 1 925t I Vitensel 1 95 31) f or de scribing convec- 

tion w ithin the atmosphere. For our spectra, we used the formula provided bv lLudwig. Frevtag & Steffen 
(119991) . which has been calibrated using 3D RHD models: 



a-aQ + {ax+ (03 + 05^^ + atgs)Ts + a^gs)Ts + a2gs 



(2) 



with 



reff-5 770K liQiog(9) 

Ts = t^ttt:: and gs = log 



1000 



and coefficients given by: 



27 500 



fl() 


fli 


fl2 


fl3 


fl4 


fl5 


ae 


1.587 


-0.054 


0.045 


-0.039 


0.176 


-0.067 


0.107 



(3) 



For matching synthetic spectra with observed ones, we need micro-turbulence as an addi- 
tional adhoc parameter Our definition of this is that a large scale (macro-) turbulent motion 



A New Library of Synthetic Spectra from PHOENIX and its Application to Fitting MUSE Spectra 















Micro turbulence [km/s] 








6 




■ 








■ 










^ 1 





ISiiiSI 








!b:1 : 




■^■■■^^■1 






■^■■■^■i 




1 









Figure 2. Distribution of micro-turbulences for that part of the grid with solar abundances for different 
effective temperatures Teff and surface gavities log(g). Color-coded is the micro-turbulence from Okm/s 
(black) to 6km/s (white). 



triggers a small scale (micro-) turbulent motion on len gth scales b elow the photon main free 
path length, which affects the strength of spectral lines (IGravl 120051) . In this picture, the micro- 
turbulence is strongly related to macro-turbulent motion, therefore we use u„iicro - 0.5 • (Vcom) 
as an ex perimental formula that follows fr om 3D radiative hydrodynamic investigations of cool 
M-stars dWende. Reiners & Ludwidl2009 ). So we first calculate the model atmosphere and then 
synthesize a spectrum from it using a micro-turbulence, which is assumed to be half the mean 
convective velocity in the photosphere. Fig.|2]shows the distribution of micro-turbulences in our 
grid for solar abundances. Unfortunately we had a problem with PHOENIX concerning convec- 
tion for giants around 7 000 K, so we had to disable convection for those models. Therefore there 
is no micro-turbulence included as well. 



4. Fitting stellar parameters 



The spectroscopic analysis of crowded stellar fields, such as star clusters or nearby galaxies has 
been limited to relatively small samples of stars thus far. The main problem in this respect is that 
traditionally used techniques like multi-object spectroscopy are restricted to the blighter, isolated 
stars in the field. We have developed a new method to overcome this limitation using integral 
field spectroscopy. Taking advantage of the combined spatial and spectral coverage provided 
by an integral field s pectrograph, we developed a riew analysis appro ach which we call "crowded 
field spectroscopy" f Becker. Fabrika & Rothl2004t iKamannlin prep .*). Via PSF fitting techniques, 
single object spectra for the stars above the confusion hmit are extracted. This deblending tech- 
nique works so well that we obtain clean stellar spectra for a significantly higher number of stars 
than hitherto possible. 

For the extracted spectra we determine the stellar parameters using a weighted constrained 
non-linear least-squ ares minimization (Levenberg-Marqu ardt), simi l ar to t he ULySS package by 
Koleva et alj ( 120091) . which has been used as well e. g. bv lWu et al.l (l201lb . Usually six different 
parameters are fitted: effective temperature T^ff, surface gravity log(gr), metalhcity [Fe/H], a- 
element abundance [a/Fe], radial velocity Uiad and line broadening cr. 

In every iteration of the Levenberg-Marquardt algorithm, a model spectrum is extracted from 
the PHOENIX grid using an N-dimensional spline interpolator. Then a line of sight velocity 
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Figure 3. The color-magnitude diagram for a simulated MUSE data cube is shown on the left. For the 
analysis we only used 1/9 of the cube, for which the stars are marked in red. On the right the errors in fitted 
radial velocities are plotted. 



distribution (LOSVD) is applied to the model, which adds line broadening cr and a shift caused 
by e. g. radial velocity v^^^. Finally an N-dimensional Legendre polynomial (i. e. the continuum 
difference between model and observation) is determined using a linear fit in a way that the model 
multiplied with this polynomial matches the observation as best as possible. Due to this we are 
independent of the continuum and the fit is done on spectral lines only. 

Figure |3] shows some first preliminary res ults for a simulated M USE data cube based on 
real HST observations of 47 Tuc, obtained by ISaraiedini et al.l (l2007h . As one can see, at the 
main-sequence turnoff point at ~17mag we still can fit the radial velocity with an accuracy of 
about 5km/s. The systematic offset in the results is caused by a known error in the creation of the 
simulated data cube. At that magnitude, the corresponding error in the fitted effective temperature 
is of the order of 100 K. 



Kirbv. Guhathakurta & SnedenI (12008 ^ showed that even with medium resolution spectra it is 



possible to determine the abundance of single elements, in their case it was some of the alpha 
elements (Mg, Si, Ca, Ti). For this analysis, they created a mask in order to fit only those parts of 
a spectiTim against a grid of synthetic spectra, where it changes most when varying the analysed 
element. They created the mask by looking for regions, where two spectra with a given Teff, a 
fixed log(g'), a metallicity of {FefH} = -1.5 and element abundances of {XIFe\ - +0.3 differ 
by more then 0.5%. Masks for several different temperatures where combined into a single mask 
that was used for the analysis. The main advantage of this method was that they did not need new 
dimensions in the grid for every new element, but could use an existing one (here {ajFe}), since 
the masks did not overlap. 

With our new PHOENIX grid, we can go one step further and create a mask specifically for 
ev ery single sp ectrum that we want to analyse. Therefore we use the same method as described 
bv lKirbv et alj (12008.) . but use the previously fitted values for Teff and \og{g) of the observed spec- 
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trum. Using the individual mask for each star, alpha element abundances can then be determined 
with an uncertainty of typically 0. 1 -0.2dex. Of course our intention is to extent this method even 
further in order to fit the abundance of other elements. 

When observing the same field multiple times, we will have radial velocities for several 
epochs, so that we can determine orbital parameters for binaries that we find. Furthermore we 
can extend our method described above to fit simultaneously the two components of a binary star 
and henceforth derive the atmospheric parameters of both stars. 

5. Conclusion 

We presented a new extensive grid of synthetic stellar spectra from PHOENIX atmospheres with 
a wavelength range and resolution that should cover all existing and upcoming instruments. Cur- 
rently its parameter range is optimized for the analysis of globular clusters, but we intend to 
extend it to higher temperatures. We also introduced a new self-consistent way of describing the 
micro-turbulence in model atmospheres. 

Furthermore we presented a first view on the methods for analyzing globular clusters with 
data obtained with the VLT MUSE together with preliminary results on simulated data. We 
showed that we will be able to examine both the kinematics as well as the binary fraction. In 
addition we will have accurate stellar parameters for most of the stars in the field. 
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